Large Scale Quantum Computation in an Anharmonic Linear Ion Trap 
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We propose a large-scale quantum computer architecture by stabilizing a single large linear ion 
chain in a very simple trap geometry. By confining ions in an anharmonic linear trap with nearly 
uniform spacing between ions, we show that high-fidelity quantum gates can be realized in large lin- 
ear ion crystals under the Doppler temperature based on coupling to a near- continuum of transverse 
motional modes with simple shaped laser pulses. 
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Trapped atomic ions remain one of the most attractive 
candidates for the realization of a quantum computer, 
owing to their long-lived internal qubit coherence and 
strong laser-mediated Coulomb interaction []], 0, H, H|. 
Various quantum gate protocols have been proposed 
0, [H, @, 0, B 0] and many have been demonstrated with 
small numbers of ions 0, H3, [n|, E3, [H, Q. The cen- 
tral challenge now is to scale up the number of trapped 
ion qubits to a level where the quantum behavior of 
the system cannot be efficiently modeled through clas- 
sical means [4|]. The linear rf (Paul) trap has been the 
workhorse for ion trap quantum computing, with atomic 
ions laser-cooled and confined in ID crystals [H, Q, 0, Q 
(although there are proposals for the use of 2D crystals 
in a Penning trap [15] or array of microtraps [6[). How- 
ever, scaling the linear ion trap to interesting numbers of 
ions poses significant difficulties 0, 0] • As more ions are 
added to a harmonic axial potential, a structural insta- 
bility causes the linear chain to buckle near the middle 
into a zigzag shape [l6[ , and the resulting low- frequency 
transverse modes and the off-axis rf micromotion of the 
ions makes gate operation unreliable and noisy. Even in 
a linear chain, the complex motional mode spectrum of 
many ions makes it difficult to resolve individual modes 
for quantum gate operations, and to sufficiently laser cool 
many low-frequency modes. One promising approach is 
to operate with small linear ion chains and multiplex the 
system by shuttling ions between multiple chains through 
a maze of trapping zones, but this requires complicated 
electrode structures and exquisite control of ion trajec- 
tories [lEzt 

In this paper, we propose a new approach to ion quan- 
tum computation in a large linear architecture that cir- 
cumvents the above difficulties. This scheme is based on 
combination of several ideas. First, an anharmonic axial 
trap provided by static electrode potentials can stabi- 
lize a single linear crystal containing a large number of 
ions. Second, tightly-confined and closely-spaced trans- 
verse phonon modes can mediate quantum gate opera- 
tions in a large architecture [18j, while eliminating the 



need for single-mode resolution and multimode sideband 
cooling. Third, gate operations on the large ion array 
exploit the local character of the laser-induced dipole in- 
teraction that is dominated by nearby ions only. As a 
result, the complexity of the quantum gate does not in- 
crease with the size of the system. 

The proposed ion architecture is illustrated in Fig. dJ 
It is a large linear array where the strong confinement 
in the transverse (x, y) direction is provided by the pon- 
deromotive Paul trap with an effective potential of the 
form V(x,y) = (mco^x 2 + muOyy 2 ^ /2, where m is the 
mass of each ion. The ions are initially Doppler cooled, 
with a number of ions at the edges of the chain contin- 
uously Doppler cooled in order to overwhelm any heat- 
ing that occurs during the gate operation. The middle 
portion and majority of the ion chain is used for quan- 
tum computation. Given an appropriate axial static po- 
tential V(z) from the axially-segmented electrodes, we 
assume these computational ions are distributed nearly 
uniformly, with a neighboring distance of about ~ 10 
jim. This enables efficient spatial addressing with fo- 
cused laser beams along the transverse direction for quan- 
tum gate operations. 

When the axial potential takes the conventional har- 
monic form V(z) = mcc^z 2 /2, the ion array is sub- 
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FIG. 1: (Color online) (a) Linear architecture for large scale 
quantum computation, where lasers address individual ions 
and couple to local modes of the ions, while the edge ions are 
continuously Doppler laser-cooled. 
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ject to the well-known zigzag transition unless the trap 
anisotropy is at least w x ,y/u z > 0.77N/^/\ogN 0, Qjf, 
where N is the number of ions. This structural instabil- 
ity is caused by the spatial inhomogeneousity of the ion 
distribution. When the ions are instead uniformly spaced 
by neighboring distance do, it is easy to see that the lin- 
ear structure is always stable even for an infinite chain 
so long as co^ y > 7C(3)e 2 /(2ra<ijj) « 4.2e 2 /(radjj), where 
£(Z) is the Reimann Zeta function and e is the charge of 
each ion. Therefore, a large linear structure can be sta- 
bilized so long as appropriate static potentials from the 
trap electrodes to make the ion distribution homogeneous 
along the axis. 

To illustrate the general method, here we consider 
an explicit example with a quartic potential V(z) = 
ot2Z 2 /2 + a^z A /4 that can be realized with a simple five- 
segment electrode geometry as shown in Fig. [2^. Under 
a quartic trap V(z), the axial equilibrium position Z{ of 
the ith ion can be obtained by solving the force balance 
equations dU/dzi = 0, where U = ^ + V( x i V)] + 

E i<? -e 2 /k, 



Tj | is the overall potential including the 



ions' mutual interactions. We optimize the dimension- 
less ratio B = \a2/e 2 \ 2 / 3 (o^/o^) characterizing the axial 
potential to produce a nearly-uniformly spaced crystal. 
To be concrete, we consider an array of 120 ions, with 
10 ions at each edge continuously laser cooled and 100 
qubit ions in the middle for coherent quantum gate oper- 
ation. We solve the equilibrium positions of all the ions 
under V(z) (see appendix A) and minimize the variance 

for the 



in ion spacing s z = yj ^ J^=ii ( Az n ~ ^z n 
qubit ions, where Az n is the distance between the nth 
and (n + l)th ion in the chain and Az n denotes its av- 
erage. The variance in spacing is shown in Fig. [2)3 as 
a function of the parameter B. The value of s z is fairly 
insensitive to B and reaches a minimum when B « —6.1. 
Here, the distribution of ion spacings z n is shown in Fig. 
Et, which is remarkably homogeneous for the qubit ions 
even though we have optimized just one control param- 
eter: s z /Az n deviates by only 3% over the entire crys- 
tal. In this configuration, if we take Az n = 10 fim for 
atomic Yb + ions, we only need a transverse center-of- 
mass frequency u Xjy /27r > 221 kHz to stabilize the linear 
structure. In this paper, we actually take uo x /2tt = 5 
MHz, as is typical in experiments, and such transverse 
confinement would be able to stabilize linear chains with 
thousands of ions under an optimized quartic potential. 

We now describe quantum gate operations with this 
large ion chain, mediated by many transverse phonon 
modes. Given the equilibrium positions of the ions, we 
can efficiently determine all axial and transverse phonon 
modes. We then apply a spin-dependent laser force, with 
the resulting interaction Hamiltonian [2(J Hl[ 

H = ^2 h^n(t)(j z n cos(kq n + fit), (1) 

n 

where the transverse displacement q n of the nth ion in the 



x direction is expressed in terms of phonon modes a& with 
eigenfrequency u k and the normal mode matrix b k by 
Qn = E/c b k n ^/h/2muo k (a\e iuJkt + a k e~ iuJkt ). The normal 
mode matrix 6^ and its eigenfrequency uu are determined 
by solving the eigen-equations E n ^ n ^n = ^V 3 ^ where 
Ai n = d 2 U/dxidx n are calculated at the ions' equilib- 
rium positions Z{. In Eq. (1), o z n is the Pauli spin opera- 
tor for the nth ion, Q n (t) denotes the Rabi frequency of 
the laser pulse on the nth ion with detuning \i from the 
qubit resonance, and the effective laser momentum kick 
k is assumed to be along the transverse x direction. (For 
twin-beam stimulated Raman laser forces and hyperfine 
state qubit s, the effective laser kick is simply along the 
difference wavevector ki — k 2 of the two beams.) Due to 
the strong transverse confinement, the Lamb-Dicke pa- 
rameter r]k = \k\ y/h/2mujk <^ 1, and the Hamiltonian 
H can be expanded as H = — E n k ^Xn(t)g k (a\e lUJkt + 
a k e~ iuJkt )a^ with g k n = r] k b k n and Xn(f) = Q n (t) sin(/it) 
(the effect of higher-order terms in the Lamb-Dicke ex- 
pansion will be estimated later). The corresponding evo- 
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FIG. 2: (Color online) (a) Sample five-segment linear ion trap 
with voltages V%(i — 1,2) to produce a quartic axial poten- 
tial. The ions are confined in the central segment, (b) The 
variance of the ion spacings s z in a linear quartic trap as a 
function of the trap parameter B that characterizes the ratio 
of quadratic to positive-quartic nature of the potential, (c) 
The distribution of the ion spacings at the optimum value 
B = —6.1. The computational ions are within the dashed 
lines, where the spacing is essentially uniform. 
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lution operator is given by [2l| 



U (r) = exp 



n.k 



[^(t)4+^*(tkk 



(2) 



where «n( r ) = Jo Xn(t)9n etUJkt ^ charac- 
terizes the residual entanglement between 
ion n and phonon mode k and (j) m n(^) = 
2 Jq dt 2 /J 2 dh Y, k 9m9nXm(t 2 )Xn(ti) sin u; k (t 2 - ti) 
represents the effective qubit-qubit interaction between 
ions m and n. 

For a two-qubit gate on an ion pair i and j, we di- 
rect laser light exclusively on these two ions (Qi(t) = 
Qj(t) = Q(t) and all other Q n (t) = 0), and the evolu- 
tion operator reduces to the standard controlled 7r-phase 
(CP) gate for a^(r) = and (j)j n {r) = 7r/4. For a large 
ion crystal, the residual entanglement with the motional 
modes cannot be eliminated completely, but we can min- 
imize the resulting gate infidelity by optimizing the laser 
pulse shape Q(t) [21]. Assuming each phonon mode k 
is cooled to temperature the infidelity of the CP 
gate from the residual motional entanglement is given 
by 5F= [6 -2(^ + ^0 -r+-r_]/8 [21], where T i{j) = 
exp[-E fe l<-)(r)| 2 A/2], r± = exp[-EJ^(r) ± 
a)[r)\ 2 p k l% and p k = coth(hw k /k B T k ). 

To minimize the gate infidelity 5F, we break the laser 
pulse on the two ions into uniform segments of constant 
intensity as shown in Fig. [3)3 and optimize the values 
(i = 1, . . . , M) over M equal-time segments [2l|. The 
control of Q (t) from one segment to the next can easily 
be accomplished with optical modulators. (Alternatively, 
we can modulate the detuning \i of the laser pulse as a 
control parameter.) After optimization of fiW, the infi- 
delity SF is shown as a function of the detuning \i in Fig. 
[3^ for M = 5 segments. For this example, we perform 
a CP gate on the 59th and 62 th ions in this 120-ion 
chain. With an appropriate choice of fiW and /i, the 
infidelity can be made negligible (well below 10 -5 ). For 
this calculation, we assume Doppler cooling for all modes 
and take gate times r in the range 50to to 500to where 
ro = 2tt/uj x = 0.2/xs is the period of transverse harmonic 
motion. The gate can certainly be faster with stronger 
laser beams (there is no speed limit), and with a faster 
gate, the control becomes easier as the gate becomes more 
localized (Fig. EJ. 

Interestingly, we use only a few control parameters 
(M = 5 segments) to perform a high-fidelity gate that 
involves excitation of hundreds of transverse phonon nor- 
mal modes. This is possible because the gate has a local 
character where the contribution to the CP gate comes 
primarily from the spin-dependent oscillations of the ions 
that are close to the target ions. To show this, we plot 



the response of each ion in Fig. [3k during the gate op- 
eration. Note that the displacement q n of the nth ion 
is spin-dependent during the gate, and we can use its 

largest magnitude over the gate time r to charac- 

terize the response of ion n, as is shown in Fig. [3fc. The 
ion response decays very fast from the target ions (59th 
and 62th in this case) and can be safely neglected after a 
distance of a few ions. Thus during a gate, only the mo- 
tion of ions near the target ions are important, and the 
other ions largely remain in their equilibrium positions. 
The resultant control parameters from this approxima- 
tion are almost identical to those shown in Fig. [3)3. Ow- 
ing to the local character of the gate, the complexity of a 
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FIG. 3: (Color online) (a) The gate infidelity SF as a func- 
tion of the the laser detuning /x from the qubit resonance, 
with optimized Rabi frequencies over M = 5 equal segments 
of the laser pulse. Different curves correspond to different gate 
times. The 120 transverse phonon modes are all distributed 
within the narrow frequency range indicated by the two ver- 
tical lines, (b) The shape of the laser pulse that achieves the 
optimal fidelity (with £F = 8.5xlO -6 ),at the detuning shown 
by the arrow in Fig. 3a (with (fi — uj x )/uj x = 9.3 x 10 -3 ) and 
the gate time r = 500to. The dashed (or dotted) lines rep- 
resent the approximate optimal solutions of the laser shape 
where only 4 (or 8) ions (from the 59th to 62th or the 57th 
to 64th, respectively) are allowed to vibrate and all the other 
ions are fixed in their equilibrium positions. This approxi- 
mation doesn't significantly change the optimal laser pulse 
shape compared to that of the exact solution (represented by 
the solid line) where all ions are allowed to vibrate, so the 
gate is essentially local and the gate complexity does not de- 
pend on the crystal size, (c) The relative response of the ions 
for the gate shown in Fig. [3]3 (characterized by the largest 

spin-dependent shift ^q^^ during the gate time r). We take 

a relative unit where | q^ | for the target ions haven been 

normalized to 1. The fast decay of the response as one moves 
away from the target ions (59th ad 62th) shows that the gate 
involves vibration of only local ions. 
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gate operation does not depend on the chain size, and we 
can perform gates in parallel on ions in different regions 
of a large chain. 

We now discuss several sources of noise for gates in a 
large ion crystal and show that their effects are negligible. 
First, the axial ion modes have large phonon occupation 
numbers under Doppler cooling alone, and the result- 
ing thermal spread in position along the axial direction 
can degrade the effective laser interaction. For exam- 
ple, the lowest axial mode in a 120-ion chain of Yb + 
ions with a spacing ~z ~ 10/im has a frequency of only 
colo/2tt = 9.8 kHz and a mean thermal phonon number 
no ~ ~ 10 3 under Doppler laser cooling (radiative 

half-linewidth j/2tt = 10 MHz). We assume the quan- 
tum gate laser beams are directed along the transverse 
direction with an axial Gaussian laser profile Q (z) oc 
e -(z/w) centered on each ion. The beam waist is taken 
as w = z/2.5 « 4/im so that the cross-talk error prob- 



-2(z/w) 2 



< 10" 



ability between adjacent ions is P c - 
The position fluctuation Sz n of the nth ion causes the 
effective Rabi frequency to fluctuate, resulting in a gate 
infidelity 5F 1 « (tt 2 /4) (5Sl n /Cl n ) 2 « (tt 2 /4) (Sz n /w)\ 
The fluctuation 5z n can be calculated exactly from sum- 
mation of contributions of all the axial modes and its 
value is almost independent of the index n for the com- 
putational ions (see Appendix). Under Doppler laser- 
cooling, Sz n « 0.26/im and the corresponding infidelity 
is SFi = 4.4 x 10 -5 . The position fluctuation of the 
ions may also lead to anharmonic ion motion, whose 
contribution to the gate infidelity can be estimated by 
SF 2 ~ (Sz n /z) 2 ~ 6.8 x 10 -4 . Finally, in the trans- 
verse direction we estimate the infidelity caused by the 
higher-order expansions in the Lamb-Dicke parameter. 
As all the transverse modes have roughly the same fre- 
quency uok ~ oj x , the effective Lamb-Dicke parameter for 
the transverse modes is r] x | Ak| ^h/2muj x « 0.038 for 
Yb + ions at lu x /2tt = 5 MHz, with each mode containing 
a mean thermal phonon number n x « 2.0 under Doppler 
cooling. The resultant gate infidelity is estimated to be 
SF 3 w n 2 4 (n 2 + n x + 1/8) w 7 x 10~ 4 [3, Note 
finally that sideband cooling is possible in the transverse 
direction as all the modes have nearly the same frequency, 
thus reducing the gate infidelity due to transverse ther- 
mal motion by another order of magnitude. 

In summary, we have shown through explicit examples 
and calculations that it is feasible to stabilize large linear 
ion crystals where the gate complexity does not increase 
with the size of the crystal and the gate infidelity from 
thermal fluctuations can be made negligibly small under 
routine Doppler cooling. The results suggest a realistic 
prospect for realization of large scale quantum computa- 
tion in a simple linear ion architecture. 

Appendix: Thermal position fluctuation of the 
ions along the axial direction 

We address individual ions through focused laser 



beams which typically take a Gaussian shape along the 
z direction with Q n (z) oc e~ z ™ / A , where z' n = z — z n 
is centered at the equilibrium position z n of the nth ion. 
Under the Doppler temperature, the ions have significant 
thermal fluctuation of their positions along the z direc- 
tion, which leads to an effectively fluctuating laser ampli- 
tude Q n (z) and induces infidelity to the gate operation. 
This position fluctuation influences both the single-bit 
and the two-bit operations in the same way. To quantify 
the gate error caused by this fluctuation, let us consider 
a spin-flip gate operated on the nth ion as a typical ex- 
ample. For a spin-flip with a tt— pulse, the gate fidelity is 
given by F t = sm 2 (Q n r + 5Q n r) * 1 - (ir 2 /4)(5Q n /Q n ) 2 , 
where Q n is the expectation value of the Rabi frequency 
which satisfies Cl n r = tt/2 for a spin- flip gate, and 5Q n is 
its fluctuation caused by the position fluctuation of the 
ion. From Q n (z) oc e~ z ™ / A « 1 — z /2 /A 2 around the 
equilibrium position, the gate infidelity SFi = 1 — F\ = 

, 2 \ 1/4 

(7r 2 /4) (Szn/A) 4 , where 5z n = (z^ - z' 2 ) character- 
izes the thermal position fluctuation of the nth ion along 
the axial direction. 

The phonon modes are in thermal equilibrium un- 
der the Doppler temperature T , with their den- 
sity operator given by p m = H k Y,{n k } P k \ n k) (nk\, 



where Pi- = 



V (n>k + l) nfc+1 is the probability of 
having n& phonons in the kth mode, and n& = 
ksT '/ ' (tvjOk) is the average phonon number. From z' 2 = 
tr(^ 2 p m ) and z'£ = tr(z^ 4 p m ), we explicitly have 



5z. 



1/2 



where 



yiEk{K' k ) (2n fe + l)M, fe j 
uo z ^k and b^ k denote th eigen-frequencies and eigen- 
matrices of the axial modes. For our example of a 120-ion 
chain with the Doppler temperature hsT/h = 62MHz 
for the Yb+ ions, the resultant position fluctuation Sz n 
is plotted in Fig. [4] for all the ions. One can see that for 
the computational ions (n from 11 to 110), Sz n ~ 0.26 
fim with its value almost independent of the ion index. 
The position fluctuation is still significantly smaller than 
the ion spacing (« 10 /ira), which ensures a tiny gate 
infidelity 5F\ as discussed in the main text. 
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FIG. 4: The axial position fluctuation 5z n is plotted along 
the ion chain, which is about 0.26/im for the computational 
ions (n from 11 to 110). 



